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ABSTRACT 

Most stars are born in clusters and the resulting gravitational interactions between cluster members may 
significantly affect the evolution of circumstellar discs and therefore the formation of planets and brown dwarfs. 
Recent findings suggest that tidal perturbations of typical circumstellar discs due to close encounters may 
inhibit rather than trigger disc fragmentation and so would seem to rule out planet formation by external tidal 
stimuli. However, the disc models in these calculations were restricted to disc radii of 40 AU and disc masses 
below 0. 1 M© . Here we show that even modest encounters can trigger fragmentation around 100 AU in the sorts 
of massive (~ 0.5 M©), extended (> 100 AU) discs that are observed around young stars. Tidal perturbation 
alone can do this, no disc-disc collision is required. We also show that very-low-mass binary systems can form 
through the interaction of objects in the disc. In our computations, otherwise non-fragmenting massive discs, 
once perturbed, fragment into several objects between about 0.01 and 0.1 M©, i.e., over the whole brown dwarf 
mass range. Typically these orbit on highly eccentric orbits or are even ejected. While probably not suitable 
for the formation of Jupiter- or Neptune-type planets, our scenario provides a possible formation mechanism 
for brown dwarfs and very massive planets which, interestingly, leads to a mass distribution consistent with the 
canonical substellar IMF. As a minor outcome, a possible explanation for the origin of misaligned extrasolar 
planetary systems is discussed. 

Subject headings: binaries: general — open clusters and associations: general — stars: formation — stars: 
low-mass, brown dwarfs — planetary systems: protoplanetary disks — 



1. INTRODUCTION 

The origin of brown dwarfs (BDs, to which we also include 
very low-mass stars < O.IM©) as well as the formation of 
massive exoplanets is contentious. The most common idea is 
that BDs form like stars, and planets form by accretion onto a 
core formed from dust conglomeration within a circumstellar 
disc. However, there is increasing evidence that an alterna- 
tive formation scenario for BDs and, possibly, some of the 
most massive exoplanets is required. BDs differ from low- 
mass stars in several ways. Firstly, there is the 'brown dwarf 
desert' - a lack of BD companions to stars, especially at sep- 
arations below about 30 AU within which companion planets 
and stars are common. In addition, the properties of BD-BD 
binaries are very different from star- star binaries, showing a 
significant lack of systems wider than ^ 20 AU beyond which 
most stellar binaries are found. Furthermore, while typically 
every second star is a binary system only every fifth BD is a 
binary. Clearly we must ask ourselves why stars and BDs are 
so different if they form by the same mechanism. Indeed, in 
order to numerically set-up a realistic stellar and BD popula- 
tion, stars and BDs nee d to be a ccording to different algorith- 
mic pairing rules KroupaeFal] ([2003 ). In other words, they 
follow different formation channels. 

The fragmentation of extended circumstellar discs is one 
of the most promising alternatives to th e standard scenari o 
of star-like formation for BDs (Goodwin & Whitworthll2 007'). 
Research on circumstellar disc evolution has made great 
progress in recent years, but there are still ni a ny unan- 
swered questions (iPapaloizou & TerquemI 120061: iHenningI 
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120081: iHillenbra nd 2008\ One issue that has gained sur- 
prisingly little attention in the past is the role of gravita- 
tional interactions between the disc and passing stars within 
the young host stellar cluster. Since stars are generally born 
in clusters (Forgan & Rice 2009), rather than as isolated ob- 
jects, star-disc and d i sc-disc interac tions (Boffin etal. 1998; 
IWatkins et al][l998al: IPMzner et al.l 12005: Thies et a l. 2005) 
probably play an important role for the dynamical evolution 
of such discs and thus for massive planet and BD formation 
through fragmentation of the discs. 

For a typical open star cluster hosting about 1000 stars 
within a half-mass radius of 0.5 pc, encounters closer than 
500 AU will happen approximately every 10 Myr (iThies et al.l 
2005). They are even more frequent in denser clusters in 
which are thought to form at least half of the stars in the 
Galaxy. Encounters in star clusters are therefore expected to 
play a significant role in the evolution of extended discs of 
several 100 AU radii. 

The formation of planets through core accretion is typi- 
cally expected to ha ppen within the in nermost 100 AU of 
a circumstellar disc (iHillenbrandl |2008|) . although there are 
many unanswered questions e specially concerning the very 
first stages of core accretion (Henning 2008). Fragmenta- 
tion has been proposed as an alt ernative way to fo r m the 
most massive ga s giant planets by iMayer et al.l (|2002|) : iBossI 
(ll997ll2004l ^20Q6). In recent years, however, there has been 
a growing consensus that disc fragmentation cannot produce 
gas giant planets in the inner regions of discs for two reasons: 
(1) heating from the central star stabilizes the innermost re- 
gions and thus inhibits fragmentation, and (2) these regions 
cannot c ool fast enough to allow temporary gas clumps to 
collapse (iLodato et al.l[2007l: IStamatellos & Whitworthll2008l: 



ICai et al .1 120 10*). This picture changes greatly beyond about 
100 AU where cooling becomes sufficiently efficient for 
gravitational collapse to occur. Consequently, disc frag- 



mentation has been propos ed as an alternative mechanism 
for the formation of BDs (IWhitworth & StamatellosI [20061: 
iGoodwin & Whitworthll2007|) . Most of these investigations, 
however, assume that the disc evolves in isolation. 

Star-disc and disc-disc collisions have been investigated in 
previous studies (Boffin etal. 1998; Watkins et al. 1998a.b) 
but without a realistic treatment of th e radiative heat transfer 
in the disc. From SPH computations 'Forgan & Ricd ('2009') 
deduce that close encounters inhibit fragmentation in typical 
protoplanetary discs (assuming an initial radius of 40 AU) 
rather than inducing it. In contrast to that, a recent study 
of disc-disc collisions ( Shen et al.|[2Q^1 0) has, however, shown 
that the direct compression of gas by the close encounters of 
stars hosting massive extended discs (disc mass ^ 0.5 M©, 
disc radius ^ 1000 AU), may trigger the formation of sub- 
stellar companions. Howev er such massive, exten ded discs 
are probably short-lived (Eisner & Carpenter 2006) (in con- 
trast to typical protoplanetary discs which are about ten times 
smaller and less massive), so that pairwise interactions of two 
discs of this kind are expected to be rare. The purely gravita- 
tional interaction of a low-mass star with no (or only a small) 
protoplanetary disc with a star hosting a massive disc would 
be much more likely. Although even these encounters may 
probably account for only a fraction of all BDs, especially in 
loose associations like Taurus-Auriga, they provide an addi- 
tional channel of forming BDs out of discs that would oth- 
erwise probably dissolve without ever fragmenting. In this 
paper we address this issue. In Section [21 besides an estimate 
of the encounter probability, the model of the circumstellar 
disc is described. Section [3| depicts the computational meth- 
ods and gives an overview of the model parameters studied 
in this work. The results are then presented in Section [4] and 
discussed in Section [5] 

2. MODEL BASICS 
2.1. Encounter probability 

The likelihood of a close stellar encounter with an impact 
parameter b depends on the stellar density of the star- forming 
environment and the impact parameter itself. The encounter 
event is then further characterized by the relative velocity of 
the two stars and thus the eccentricity of the path. The conse- 
quences on the circumstellar disc further depend on the mass 
of the perturber and the host star (i.e., the ratio of binding en- 
ergy and magnitude of the perturbation) and the inclination of 
the disc plane wrt. the plane of the encounter hyperbola. 

Stars can be born in a variety of environments from low- 
density groups like Taurus-Auriga (about 300 stars distributed 
in groups, each being <1 pc in radius and containing about 
a few dozen stars) up to ONC-type clusters (thousands of 
stars within 0.5 pc) or even in more massive and dense clus- 
ters (Kroupa 2005). An estimate of the encounter probability 
has been derived assuming a Plummer star cluster model of 
half-mass radius, ri/2, total mass Mtot, and an average stel- 
lar mass, M^. Note that two other radial scales are also of- 
ten used to describe a Plummer model: the Plummer radius, 
rpi = V2^/^-lri/2 ^ 0.766ri/2, and the gr avitational s cale ra- 
dius, Tgrav ^ 2.602ri/2 ^ 3.395rpi. From Thies et al.l ([2005), 
Section 1.3, the expected time, ^enc, between two encounters 
within an encounter parameter, /?, can be obtained from the 
characteristic crossing time ^cr, 
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Fig. 1. — The Toomre parameter Q as a function of the radius from the 
central star for both the initial disc setup (solid line) and the settled disc after 
20 kyr (dashed line). The profiles of both disc stages are nearly identical 
below about 200 AU while Q is slightly larger (i.e. the disc is more stable) 
between 200 AU and 400 AU, the outer rim of the initial disc. The region 
within 50 AU is skipped since it is initially partially void. 



and the average number, ^enc, of encounter parameters below 
b, for a given number of stars, N, 
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It has to be noted that the actual periastron distance, rperi, is 
always (and sometimes significantly) smaller than b due to 
gravitational focusing. The relation between rperi and b for a 
given eccentricity, e, is 
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For an ONC-type cluster with N = 7500, an average star mass 
M^ = O.5M0, a total mass Mtot = 10000 M© (i.e 30%ofthe 
mass consists of stars and the rest of gas, see iKroupa et al.l 
|2001|) . and a half-mass radius of ri/2 = 0.5 pc we find that en- 
counters below 500 AU occur on average every 11 million 
years. If a massive extended disc exists for about 1 Myr this 
means that about 9 per cent of all d iscs of this type will suffer 
such an encounter. See also lThies~ et al. ( 2005i) . 

2.2. Disc model 

The initial conditions for the di sc are taken from 
IStamatellos & WhitworthI (I2008L 12009a ). The disc model has 
a po wer law profile for temperature, T , and surface density, S 
from IStamatellos & WhitworthI (l2008h . both as a function of 
the distance, R, from the central star. 
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where q^ = 1.75, qr = 0.5, and Tq = 300 K. So is chosen cor- 
responding to a disc mass of 0.48 and 0.50 Mq (see Table [B, 
i.e., So = 8.67 x 10^ gcmT^ and So = 9.03 x 10^ gcmT^, re- 
spectively. Too = 10 K is a background temperature to account 
for the background radiation from other stars within the host 
cluster. This co nfiguration lea ds to an initial Toomre stability 
parameter, Q, (lToomrdri964)) between 1.3 and 1.4 between 
about 100 and 400 AU, leading to weak spiral density patterns 
but not to fragmentation. Q is given by 
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with Cs being the sound speed, and Hi being the epicyclic fre- 
quency which can be, at least roughly, approximated by the 
Keplerian orbital frequency. The radiation of the perturber 
is included in the same way as that of the central star, as- 
suming a mass-to-luminosity relation of L oc M^ for low-mass 
stars. However, as long as the perturber is no more than about 
0.5 Mq, corresponding to an early M-type or late K-type main 
sequence star, its influence on disc evolution is very small. 

The disc is initially populated by SPH particles between 
40 and 400 AU. Before starting the actual encounter compu- 
tation the disc is allowed to "settle" for about 20 000 years 
(i.e., about two orbits at the outer rim) to avoid artefacts from 
the initial distribution function. During this settling the disc 
smears out at the inner and the outer border extending its ra- 
dial range from a few AU up to > 500 AU, eventually stabilis- 
ing. The Q value is slightly increased beyond 200 AU relative 
to the initial setting while remaining largely unchanged below 
this radius, as can be seen in Figure [T] The actual surface den- 
sity after the settling is therefore slightly smaller and thus the 
disc slightly more stable than in the initial setup. Any self- 
fragmentation of the unperturbed disc would have happened 
during the settling time when the Toomre parameter is lowest. 

3. METHODS 
3.1. Smoothed particle hydrodynamics (SPH) 

All computations were performed by usi ng the well-tested 
DRAGON SPH code by iGood win et al l (12004 includin g 
the radiative transfer extension by Stamat ellos et al.l (l2007bb . 
Most of the numerical p a rameter s ettings have been adopted 
from IStamatellos et al] (12007 ah : IStamateUos & WhitworthI 
(|2009a '). Gravitationally collapsing clumps are treated by 
sink particles which form if the volume density exceeds 
10~^gcm~^, and the clump is bound. The sink radius is 1 AU, 
in accordance with the studies mentioned above and in agree- 
ment with the local Jeans criterion. In addition, the central 
star and the perturbing star are both represented by sink parti- 
cles. 

There are ongoing discussions whether fragmentation can 
be triggered artificially due to the n umerical be haviour of 
SPH and grid-based hydrocodes. iHubber et al.l (|2006|) find 
that smoothed particle hydrodynamics ( SPH) do e s not suffer 
from artificial fragmentation. However, iNelsonI (|2006|) sug- 
gest that SPH may artificially enhance fragmentation due to 
a pressure-underestimation if the smoothing radius, h, is con- 
siderably larger than about 1/4 of a vertical scale height of 
a circumstellar disc. Stama tellos & Whitwort h (2009b) show 
that this criterion is fulfilled for 150,000 or more particles, for 
th e types of disc studied in this paper. 

ICommercon et al.l (l2008h have performed calculations of 
fragmenting protostellar clouds with the adaptive mesh refine- 
ment (AMR) code RAMSES and the SPH code DRAGON 



(which we have been using for this study) and compared 
the results at different resolutions given by the number N 
of grid cells and particles, respectively. For about A^ = 
5 X 10^ cells/particles fragmentation in both codes appeared 
nearly equal. With N = 2 x 10^ the results were very sim- 
ilar to the high-resolution AMR computations while be- 
ing still acceptable for niany p urposes at N = I x 10^ . 
IStamatellos & Whitw orth' ('2009b) briefly review the reso- 
lution criteria with respect to their models using 150 000 
in that paper and 250 000 to 400 000 particles earfier 
(IStamateUos & Whitworth..2009a) . 

The three most important resolution criteria are the resolu- 
tion of the local Jeans mass, Mj, (and thus the Jeans length), 
the local Toomre mass and the vertical scale height of the disc. 
In particular, the local Jeans mass. 
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must be resolved by at least a factor of 2 x A^neigh = 100, corre- 
sponding to 0.5 Mj for the 100,000 particle discs and 0.2 Mj 
for the 250,000 particle discs. Since the global disc setup as 
well as the evolution aft er the perturbation (Secti onB]) is quite 
similar to that used by (IStamatellos & Whitwor th 2009b), as 
similar range of Jeans masses of the forming clumps, i.e., 
~ 2-20 Mj, can be reasonably assumed for our models and 
has actually been tested for the models E/X002 and E/X009 
for all forming objects. The minimum Mj during the clump 
formation is typically between 4 and 6 Mj and does never go 
below about 2 Mj. Therefore, Mj is resolved by a safe factor 
of at least 4 in the low-resolution case and 10 in the high- 
resolution case. Accordingly, the minimum Toomre mass of 
2.5 Mj is adequately resolved as well as the vertical scale 
height (by at least a few smoothing lengths). 

3.2. Radiative transfer model 

In iThies et aP (l2005h the possibility of fragmentation has 
been estimated via the Toomre criterion for both isothermal 
and adiabatic equations of state, both resulting in perturbed 
regions with highly unstable conditions (Q < 1). However, 
Toomre instability does not necessarily lead to actual frag- 
mentation. 

In realistic models of radiative transfer, such as the one used 
here, the thermal response of the gas to density changes is be- 
tween near-isothermal in regions with efficient cooling (out- 
side about 100-200 AU) and near- adiabatic in r egions with 
long cooling times (the inner parts of the disc). Ishen et aO 
(|20 10) have shown that direct disc-disc collisions may lead 
to fragmentation of massive extended discs at some 400 AU 
even for a "thick disc approximation", i.e., the near-adiabatic 
case. 

The Stam atellos et al.l ([2007b) method uses the density and 
the gravitational potential of each SPH particle to estimate 
an optical depth for each particle through which the parti- 
cle cools and heats. The method takes into account com- 
pressional heating, viscous heating, radiative heating by the 
background, and radiative cooling. It performs well, in both 
the optically thin and opticall y thick regimes, and has been 
extensively tested by Stamate llos et al.| (|2007bl). In particular 
it reproduces the detailed 3D results of iMasu naga & Inu tsiikal 
(2000), Boss & BodenheimeB (119791) . lBoss & Myhifl (1991), 
( Whitehouse & Bate 2006\ and also the analytic results of 
■Spieg el (1957ii). Additionally the code has been tested and 



performs well in disc configu rations as it reproduces the ana- 
lytic results of lHubenvl(ll99Qh . 

The gas is assumed to be a mixture of hydrogen and helium. 
We use an equ ation of stare by Black & Bodenheimer (1975); 
iMasunagaetal. (1998); Boley et al. (2007) that accounts (1) 
for the rotational and vibrational degrees of freedom of molec- 
ular hydrogen, and (2) for the different chemical states of hy- 
drogen and helium. We assume that ortho- and para-hydrogen 
are in equilibrium. 

For the dust and gas opacity we use the parameterisation by 
| Bell & Lin (1994), ^{p^ T) = hcq p^ T^ , where hcq, a, b are con- 
stants that depend on the species and the physical processes 
contributing to the opacity at each p and T. The opacity 
changes due to ice mantle melting, the sublimation of dust, 
molecular and H~ contributions, are all taken into account. 

3.3. Overview of Calculations 

We have conducted thirteen SPH calculations using 
100,000 and eight using 250,000 SPH particles (i.e., 21 mod- 
els in total) corresponding to about nine CPU months in total 
on 16-core machines. The model identifiers that begin with an 
"E" or "F" refer to the low resolution models while the mod- 
els with "X" or "Y" are the high resolution ones. Subsequent 
letters with identical numbers correspond to follow-up cal- 
culations (beginning at the moment of the closest encounter) 
with identical model settings. Due to the dynamical interac- 
tion of the parallel computing CPUs a slight random perturba- 
tion is imposed to the subsequent evolution such that the out- 
come differs within the statistical noise. The high resolution 
calculations are set up with parameters of preceding lower- 
resolution models that showed fragmentation. Therefore, the 
fact that all high resolution computations show fragmentation 
is due to the parameter selection. The low-resolution calcula- 
tions are used as a parameter survey while the high resolution 
ones are follow-ups to selected low-resolution ones. Please 
note that not all non-fragmenting models are listed here but 
only those which represent the borders of the parameter space 
between fragmentation and no fragmentation. 

We chose a 0.5 M(7) p erturbing star as a typical member of a 
star cluster (Kroupa" 2001) for the majority of models, but we 
also performed calculations with perturbers down to 0.3 Mq. 
The encounter orbit is slightly inclined (^ = 10° against the ini- 
tial disc plane) with varying eccentricities from e=lA (near- 
parabolic) to ^ = 2 (hyperbolic), corresponding to a relative 
velocity of 0.4-2.7 km s~^ at infinity. For these parameters the 
disc typically fragments to form a few very-low-mass stars 
and substellar-mass objects. A calculation with e = ?> has been 
performed but yielded no fragmentation. The same holds true 
for encounters at 600 AU or more. Depending on the eccen- 
tricity the calculations start 5000-10000 years before perias- 
tron to ensure a sufficiently low initial interaction. The calcu- 
lations continue for 15 000 years after encounter for the mod- 
els E/X/Y002 and at least 20 000 for the others. The masses 
and orbital radii of the objects formed in the calculation are 
determined at this time. After 15 000 years the accretion pro- 
cess of the objects has largely finished while dynamical evo- 
lution may still lead to major changes of the orbital param- 
eters (and even to ejections of some bodies) if the calcula- 
tions would have been continued over a longer time interval. 
For this reason, the orbital separations computed in this study 
have to be taken as a preliminary state. 

It has to be noted that the orbit is altered slightly during 
the passage due to the transfer of mass and angular momen- 
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Fig. 2. — Snapshot of a circumstellar disc modelled with 250,000 SPH par- 
ticles, around a Sun-type star being perturbed by a close star-star encounter 
(model X002, Table [T). The time stamp in each frame refers to the time of 
the encounter. 



tum. Furthermore, since the encounter dissipates energy the 
post-encounter speed is typically lower than the pre-encounter 
speed and in some cases both stars may be captured in an ec- 
centric binary (which actually happened to the companions la 
and lb in calculation E003; see Table [T] and Figures [3] and lU). 
We found, however, that these effects are small. 

4. RESULTS 

4.1. General findings 

As the unperturbed disc is marginally stable it does not 
fragment until the passage of the perturbing star. The first 
visible effects of a typical fly-by are the appearance of tidal 
arms and a mass transfer to the passing star (which acquires 
a small disc). About 2000-3000 years after the periastron 
passage parts of the tidal arms become gravitationally unsta- 
ble, and spiral-shaped over-densities begin to form (see the 
snapshots from model X002 in Figure [3D). Some of these 
continue to contract into a runaway collapse and form bound 
objects represented by sinks (see Section O. This typically 
happens at radii between 100 and 150 AU from the central 

^ Supplementary content like movies from our calcula- 
tions can be downloaded from the AlfA download page, 
http://www.astro.uni-bonn.de/~webaiub/german/downloads.php 
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Note. — Overview of our SPH computations. The columns (from left to right) show the computation ID, the number of gas particles in thousands, the mass of 
the central star, its disc and the perturber star. Then follows the periastron distance of the encounter, the eccentricity and inclination. The four rightmost columns 
show the amount of gas captured by the perturber in a circumstellar envelope within a radius 40 AU, the number of formed companions in total and the number 
of bodies that got ejected during the calculation as well as the number of binary systems. In all models except for E015 and E016 the perturbing star has 0.5 Mq 
and passes the central star on an initially hyperbolic orbit with initial eccentricity between 1.1 and 3.0 and inclination to the disc plane between 10° and 45°. 
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Fig. 3. — Snapshots of a forming binary around about 8000 years after the 
fly-by in model X002 (see also Figure |2). The components of the remaining 
VLMS pair have masses of 0.08 Mq and 0.09 Mq, subsequently accreting 
another 0.01 Mq each. The escaping third body has 0.05 Mq and is eventu- 
ally ejected. 



Star with some clumps forming even around 200 AU. Tempo- 
rary overdensities do also occur within less than 100 AU but, 
however, dissolve quickly. This is probably due to the heat- 



ing from the central star and the less effective cooling in these 
regions as already being me n tioned in the Introduction and 
by Whi tworth & StamatellosI (l2006h : iGoodwin & WhitworthI 
(2007). At radii larger than about 200 AU, on the other hand, 
the gas density is apparently too low for gravitationally bound 
clumps for form. 

Each of the calculations typically produces between two 
and five (six on one case) low-mass objects with masses be- 
tween those of very massive planets (0.01 Mq) and very low- 
mass stars (0.13 M©). Additionally, mass is accreted onto 
the central star, and also mass is transfered to the perturber 
during the encounter. In the model shown in Figure [2l five 
objects with masses between 0.013 Mq and 0.10 Mq form; 
two of them being bound in a binary system. In addition, 
0.016 Mq of gas accretes onto the central star within 15 000 
years, and 0.0028 Mq = 2.8 Jupiter masses (Mj) of gas trans- 
fers to the perturber (into the representing sink particle) while 
about 0.012 Mq (^ 12Mj) is captured around the perturber 
as a circumstellar envelope. In our analysis, any gas that is 
present within a radius of 40 AU but outside the sink radius is 
counted as circumstellar material of the perturber and shown 
in the ninth column of Table [T] This mass capture is typical, 
ranging between 2 and 16 Mj. 

In total, 28 bodies between 0.01 and 0.15 Mq formed in 
eight high resolution calculations, while the ten 100 000 par- 
ticle calculations that showed fragmentation, yielded 30 bod- 
ies between 0.01 and 0.16 Mq (see Tables [J and [21). Thus, 
there is no statistically significant relation between the reso- 
lution and the number of formed bodies. Similarly, no clear 
trend towards lower minimum masses of the companions for 
the high resolution can be derived from the current data. The 
average mass of the lowest-mass member is 0.056 Mq for 
the 100 000 particle models while it is 0.044 Mq for the 
250 000 ones. Similarly, the average mass of the highest- 



TABLE 2 
Outcome of Individual Calculations 



Model ID 
E002 



100 000 particles 



(Mo) 



0.021 
0.080 
0.100 
0.131 



Separation 

(AU) 



332-342 
ejected 
20-30 
50-100 



Model ID 
X002 



250 000 particles 



~1~ 

2 

3 

4a 
4b 



(Mq) 



Sep. 

(AU) 



0.013 
0.046 
0.057 
0.088 
0.099 



ejected 
ejected 
155-225 
70-110 
70-110 



E003 






X003 






la 


0.013 


130-500 


1 


0.021 


50-100 


lb 


0.086 


130-500 


2 


0.022 


200 


2 


0.130 


29-42 


3 


0.055 


ejected 








4 


0.154 


130-200 


E004 






X004 






1 


0.134 


160-190 


1 


0.026 


150-220 


2 


0.162 


60-90 


2 


0.096 


50-170 








3 


1.032 


25 


E006 






X006 






1 


0.063 


35-70 


la 


0.033 


150 


2 


0.154 


120-220 


lb 


0.046 


150 








2 


0.148 


60 


E009 






X009 






1 


0.026 


16 


1 


0.049 


20 


2 


0.099 


50 


2 


0.064 


230 


3 


0.133 


140 


3 


0.073 


230 








4 


0.138 


60 


EOlO 






XOIO 






1 


0.066 


190 


1 


0.031 


190 


2 


0.117 


80 


2 


0.058 


40 


2 


0.092 


80 


2 


0.106 


40 








3 


0.128 


80 


EOll 






XOll 







1 


0.031 


110 


1 


0.073 


100-180 


2 


0.038 


300 


2 


0.107 


50-80 


3 


0.038 


700 








4 


0.048 


230 








5 


0.074 


15-20 








6 


0.098 


125 









FOll 


Y002 




2 
2 


0.090 150-400 1 0.080 

0.168 60-100 2 0.108 

3 0.137 


P^llO 

100-600 

P^50 


E015 


Additional calculations with 100 000 particles 
E016 




1 

2 
3 


0.070 ejected 1 0.093 
0.103 10-29 2 0.149 
0.134 70-101 


13-30 
55-98 



Note. — List of the companions formed during individual calculations 
within 15 000-20 000 years after the flyby, sorted by mass. The table shows 
the masses and the approximate minimum and maximum separations from 
the central star. Note that bodies 4a and 4b in X002 got bound in a VLMS bi- 
nary with a mutual distance of about 4 AU through a triple encounter while la 
and lb in E003 got bound through dissipational capture. The binary capture 
of la and lb in X006 probably involved both mechanisms. 



mass is 0.097 M© and 0.083 M©, respectively. Table E] does 
also show the separations of the objects at the end of each 
calculation. It has to be noted that these can differ largely 
from the initial separation at the moment of formation. Dy- 
namical interaction between the objects and the disc as well 
as mutual encounters push some bodies at wide and eccen- 
tric orbits or even eject them while others migrate closer to 
the star. The closest separation observed in our calculations is 
about 10 AU in model E015 after 15 000 years of evolution. 
The overall outcome is very simi lar to that of unperturbed 
self-frag menting disc as shown in IStamatellos & WhitworthI 
(l2009bh except for the lower disc mass. This is quite plau- 
sible since the main difference between the two scenarios is 
the cause for the density patterns which undergo fragmenta- 
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Fig. 4. — Snapshots of another forming binary around about 12000 years 
after the fly-by. The components of the remaining VLMS-BD pair have very 
unequal masses of 0.09 Mq and 0.013 Mq (with a mass ratio q = 0.14). In 
contrast to the binary formed in X002 (Figure [3 these bodies became bound 
due to a frictional encounter of their accretion envelopes. 



tion. While Toomre density waves are the source of fragmen- 
tation in self-fragmenting discs, the gravitationally induced 
fragmentation occurs in tidal arms. In both cases, a dense 
bar-like in the disc reaches the density limit for fragmentation 
and thus the underlying physics is essentially the same. 

Figure[5] shows the mass distribution of the n bodies created 
all calculations with 100 000 particles (top frame), 250 000 
particles (middle frame) and both combined (bottom frame). 
A power law mass function dn/dm = ^(m/M©)"" can be fit- 
ted to this distribution. This has been done in this figure for 
the substellar regime (dashed line). In the bi-logarithmic scal- 
ing, the slope of dlogn/dlogm correspon ds to 1-a due to 
the differentiation of the logarithm (see Thi es & Kroupall2007l 
for details). The linear fit to the mass function for all calcu- 
lation outputs combined appears to have a slope of 0.9;^q4 (< 
0.16 Mq) corresponding to power law index of a = -0.1;^q;4. 
If the calculations are separated by resolution, the 100 000 
particle models correspond to a = +0.1;^^-5 while the 250 000 
particle models similarly yield a = +0.1 i^Q^. The slightly flat- 
ter IMF in the separate mass functions may be interpreted 
as a weak resolution-dependence of the average created mass 
which gets smeared out in the combined mass function. How- 
ever, the difference between the slopes of 0.2 is within the 
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Fig. 5. — Mass distribution of the bodies created in the 100 000 particle 
models (top frame), the 250 000 particle models (middle frame), and both 
combined (bottom frame). The majority of companions is found in the mass 
interval between 0.08 and 0.16 Mq while most of the rest is in the substel- 
lar regime. The substellar mass function corresponds to a power law in- 
dex of a = +0.1!?;^ for the 100k models, a = +0.l![j ^ for the 250k mod- 
els, and q; = -0.1 ^4 for both combined (see text for further explanations). 
The dark and light grey-shaded regions refer to the 1 and 2 a confidence 
limits of the fit while the errorbars correspond to the 1 a Poisson errors. 
The vertical dotted line marks the hydrogen-burning mass limit of 0.075 Mq 
fChabrier & Baraffe.2000.) . 



la uncertainty The uncertainty, based on the Poisson errors 
of the log mass bins, is indicated by a shaded region. Above 
0.12 Mq there is a sharp drop in the mass distribution that 
can be tre ated as a truncation. Similar resul ts ha ve been ob- 
tained by IStamatellos & WhitworthI (I2009al) and IShen et al.l 
(120101) . Interestingly, this is also in good agree ment with 
both the sub- stellar IMF deduced by Kroupa ('2001') as well as 
with the sep arate substellar IMF deduced in Thies & Kroupa 
(l2007Ll20Q8|) . However, more calculations are needed to pro- 
vide statistically robust tests of the IMF. Furthermore, there 
is no weighting with respect to the different likelihood of the 
different encounter settings (i.e., inclination, fly-by distance 
etc.). Since the outcome of the computations does not show 
a specific dependency on the orbital parameters (except for 
whether there is fragmentation) such a weighting might in- 
troduce an artificial bias by amplifying the random noise of 
higher- weighted models. 

4.2. Binary formation 

In the same wa y as in the models of 
IStamatellos & WhitworthI (l2009a|) of isolated circum- 
stellar discs, very low-mass binary systems can also form 
in gentle three-body encounters between low-mass objects 



in the disc. Such a triple in-orbit encounter occurred during 
the X002 computation about 8,000 years after the fly-by, as 
shown in Figure [3] where a 4 AU binary system composed 
of 0.09 Mq and 0.10 Mq objects forms through a triple 
encounter. This binary remains bound to the host star in 
an orbit of about 100 AU. These separations are consistent 
with the observations of VLMS binaries according to which 
the most pro bable separation is around 3 AU (see, e.g.. 
Figure 10 in IStamateUos & WhitworthI l2009aD . Another 
binary-forming mechanism was unveiled during calculation 
E003 (see Table©. A 2 AU binary, consisting of a VLMS of 
0.09 Mq and a BD with 0.013 Mq, formed through a grazing 
encounter of circum- sub stellar accretion envelopes, which 
subsequently evolve into a circumbinary disc. The event is 
shown in a sequence plot in Figure (H In two similar events 
in run X006, two accreting clumps merged to a single one 
of about 0.03 Mq 9600 yrs after the fly-by, and 6000 yrs 
later this body and a third one got captured into a binary 
of about 5 AU separation, probably involving both triple 
encounter and grazing collision. The masses at the end of the 
calculation are 0.03 and 0.04 Mq. 

Binaries formed via these processes may later be sepa- 
rated from their host star in a subsequent encounter with an- 
other st^_as_already discussed for sin gle substell ar compan- 
ions (iKroupall 19951: lGoodwin~& WhitworthI l2007h . or, possi- 
bly, by dynamical interaction with more massive companions 
in the same system. These models produce 55 systems of BDs 
and VLMSs, and 3 binaries giving a binary of 3/55 = 0.05, 
whereby the binaries have semi-major axes of 4, 2 and 5 AU. 
Thus, while the separations are quite consistent with the ob- 
served VLMS and BD binaries, the present theoretical binary 
fraction is somewhat low. 

4.3. The role of eccentricity and inclination 

In a previous study (iThies et al.ll2005h we found that the 
strength of the perturbation increases with decreasing ec- 
centricity and inclinatio n. This agree s with results of an 
SPH parameter study bv lPfalzner et al.l (120051) . While disc- 
disc co llisions produce most objects due to shock formation 
(IShen e t al. 2010), coplanar encounters increase the effective- 
ness of the tidal perturbations which cause our discs to frag- 
ment. This is a consequence of the low relative velocity be- 
tween the perturber and the perturber- facing parts of the disc 
in coplanar encounters, and thus a longer tidal exposure time. 
While drag forces in disc-disc collisions are larger for larger 
collisional velocities the tidal force does not depend on the 
velocity, and therefore the total amount of tidally transferred 
momentum for a given flyby distance is larger for slower en- 
counters, i.e., for those with lower eccentricity. 

5. DISCUSSION 

Our computations show that the tidal perturbations of mas- 
sive (otherwise stable) discs can form massive planets and 
BDs. Previous comp utations had ass umed the interaction 
of two massive discs (IShen et al.l 12010 ). reducing the prob- 
ability of such an event dramatically. Our scenario only re- 
quires a single extended disc. This scenario can produce 
very low-mass companions at large distances from the pri- 
mary star, and may help explain recent direct detections of 
massive planets orbiting at > 100 AU, far beyond where 
core a ccretion could h ave formed t herrfl, as well as interme- 
diate (iLiu et al.1 l2002l) and distant (IStamateUos et al.1 12007 al: 

^ See also the Extrasolar Planets Encyclopaedia, |http://exoplanet.eu/ 1 



IStama tellos & WhitworthI [2 "009aj) BD companions to stars. 
Like jStamatellos & Whitwo rt^2009ab we are also able to 
form BD binaries through three body encounters within the 
disc, and in addition, through dissipative encounters. But the 
required disc mass is considerably smaller, by about 30-40 %, 
since even initially stable discs do fragment upon perturba- 
tion in our calculations. The mass function of companions 
formed by this process is in good agreement with that of self- 
fr agmentation and with th e separate substellar IMF deduced 
in iThies & Kroupal (|2007|) . However, the volume of the cur- 
rent results does not allow robust statistical tests of the nor- 
malization of this sub- stellar IMF relative to the stellar IMF. 
Nevertheless, the obtained results are in remarkable agree- 
ment both with the form of the substellar IMF and the binary 
properties in the BD and VLMS regime, although the binary 
fraction is somewhat low (Section [4. 2t . More computations 
will have to be performed at different resolutions (up to a mil- 
lion particles) to enhance quantity as well as quality of the 
data. 

5.1. How often do such encounters occur? 

Furthermore, it has to be noted that triggered fragmen- 
tation is probably responsible only for a fraction of disc- 
fragmentation events while a good fraction may be triggered 
simply by over-feeding of an accreting disc towards self- 
fragmentation. Although the borders of the parameter sub- 
space suitable for triggered fragmentation are not yet fully de- 
termined the current findings show that encounters of a higher 
inclination than 45 degrees, outside 600 AU or with eccentric- 
ities above 2 are generally unlikely to trigger fragmentation 
for the disc type studied here. The same holds true for per- 
turbers of less than 0.3 Mq. 

By combining all these limits of the parameter space one 
can estimate the fraction of random encounters with a mutual 
periastron below 500 AU that are suitable for fragmentation 
of the analysed disc type. If a characteristic velocity disper- 
sion of r^ 2kms~^ and an upper stellar mass limit of 10 Mq 
within the host cluster is assumed only about 3 % of all en- 
counters below 500 AU lead to fragmentation. However, discs 
with even only slightly larger mass or lower background tem- 
perature may be pushed much easily to fragmentation, proba- 
bly enlarging the parameter space significantly. Furthermore, 
other perturbing effects like stellar winds or supernova shock 
waves, which might influence large volumes within the host 
cluster at once, are not covered by this study, nor are the ef- 
fects of dust. 

5.2. Consequences for planet formation 

All objects formed in our calculations have masses above 
0.01 Mq or 10 Mj and are thus above the masses typically 
assumed for planets. However, we cannot rule out at the 
moment that also objects below 0.01 M© may form through 
fragmentation, especially around lower-mass host stars which 
heat the inner disc region much less than more massive ones. 
Another critical point may be the resolution limits (although 
the local Jeans mass is well resolved even in the inner disc 
region, as discussed in Section [3Jl) . 

It has to be emphasized, that planet formation probably 
typically takes place in the inner disc region through core- 
accretion (Hil lenbrand! 2008) which are less influenced by 
the perturbation unless companions formed through it migrate 
into these inner regions. There might still be significant im- 
pacts on the outcome of planet formation for stellar encoun- 
ters, though. Even temporary gravitational instabilities that 



do not collapse into a brown dwarf or planet directly may in- 
duce the formation of substellar companions down to Kuiper- 
Belt objects by induced v orticity and subsequent dust trap- 
^ping^ dBarge & Sommeriai ri995: Kl ahr & Bodenheimedl2003[ 
2006). Also the development of baroclinic vortices may 
be altered under the influence of tidal perturbations, ei- 
ther inhibiting or promoting the formation of dust aggre- 
gates. This mechanism may even work in typical protoplan- 
etary discs and is subject of our ongoing research. Accord- 
ing to solar system architecture (J Helled 11993': "EggerseTal] 
119971: F Kenvon & Bromlev 2004) a nd rad ionuc lide evidence 
(Takigawaetal. 2008; Sahiipal & Gupta 20091: iGaidos et al.1 
2009; Gounelle & M eibom 2008, however, disagree) the 
highly probable origin of our Sun in a large star-forming re- 
gion further emphasizes the importance of such scenarios. 

Another fact worth to be discussed is the capture of disc 
material by the perturbing star. As mentioned in Section |4] 
about 0.003 M© = 3 Mj are accreted by the perturber in a typ- 
ical 500 AU encounter like model E/X002, while the amount 
of accreted gas can be as large as 10 Mj, as in models E/X003 
(see Table [B- An even larger amount can be contained in a 
circumstellar disc or envelope around the perturber. Although 
being smaller than typically assumed for the Minimum-Mass 
Solar Nebular (MMSN, see Crida 2009) this may still be suf- 
ficient for the formation of Jupiter- or Saturn- type planets 
around the perturber. Since the orientation of this encounter- 
related accretion disc is not correlated to the stellar rotation 
of the perturber, this scenario may provide an explanation for 
highly inclined or even retrograde planets which have recently 
been detected (N aritaeta l. 2009; J ohnson et al. 2009). If the 
captured gas is accreted onto a pre-existing circumstellar disc 
of the perturber star this might even lead to the formation of 
planetary systems with multiple mutually inclined (or even 
retrograde) orbital planes. 

We have to emphasize at this point that gas capture from 
a massive extended circumstellar disc is only one of several 
possible gas capture scenarios, and is being observed in our 
work as a spin-off besides the main topic of this article. As 
the more general case, capture of material from any dense 
gas aggregate in the hosting star- forming region after the for- 
mation of the protostar itself may be a possible channel to 
form non-aligned planets. As the most general formulation, 
we note that the whole process of planet formation from the 
pre- stellar cloud in the context of stellar encounters in young 
star clusters to a fully established planetary system appears to 
be a discontinuous one in many cases, probably limiting the 
probability of regularly- shaped planetary systems with Solar 
system-like architecture. This issue will be investigated in de- 
tail in future work. 

One issue not treated by our calculations is the varying pro- 
tostar and disc mass during the accretion process. According 
to Machida et al. (2010) the growing disc may become tem- 
porarily unstable when the mass of the protostar is still neg- 
ligibly small ('^ 1O~^M0). In their nested-grid calculations, 
the disc fragments in the region of < 100 AU, and even within 
10 AU. However, their simulations do not include a realistic 
treatment of radiative transfer. It is subject to future investiga- 
tions whether accreting discs under perturbation may develop 
into planetary systems with Solar-type architecture. 

Furthermore, the influence of magnetic fields are com- 
pletely ignored in our calculations. iLiverts et"an ([2010!) show 
that both hydromagnetic and thermomagnetic effects may am- 
plify density waves in the disc into instability, and may pro- 



vide an effective viscosity via turbulence which may assist 
accretion onto the star. 

Observations of tidally perturbed fragmenting discs would 
surely be the best confirmation of this scenario. Due to the 
short duration of the flyby of only about 10,000 years the 
chance of such an event to be 'caught in the act' is small. 
It may, however, be possible to identify encounter-triggered 
fragmentation if observers succeed in detecting the typical 
tidal arms around the disc-hosting star as well as smaller 
amounts of non-circular filaments around a close-by star in- 
dicating it as a candidate perturber star. Such structures 
may be detectable with future high-resolution instruments like 
ALMA in star- forming regions like the ONC. Another pos- 
sibility is to target FU Orionis stars which are thought to be 
undergoing enhanced accretion, such enhanced accretion may 
be caused by the perturbation of discs by a close encounter. 

6. SUMMARY AND CONCLUSIONS 

In a series of SPH calculations we have shown that mas- 
sive (~ O,5M0), extended (> 100 AU) circumstellar discs 
can be stable when isolated, but fragment when perturbed 
by a moderately close and slightly inclined stellar encounter. 
Binaries formed in two cases via triple encounters of com- 
panions and via grazing encounters of accreting envelopes. 
This agrees with binary formation in self-fragmenting discs 
fe tamatellos & W hitworth 2009a.). We further found that the 



mass distribution of companions formed in discs is in agree- 
ment with t he canonical substellar IMF as a separate popula- 
tion (KrouD al200Tl:lKroupa & Bouvied200"3l:lThies & Kroupal 
2007, 2008). We have, however, to note that direct tidally in- 
duced fragmentation is probably only responsible for a frac- 
tion of all disc fragmentation events due to the relatively small 
orbital parameter subset that is suitable for fragmentation of 
discs of about 0.5 Mq. More massive discs are expected to 
be more prone to fragmentation even due to weak pertur- 
bations, but are also more likely to reach the limit for self- 
fragmentation. This study, however, shows that tidal pertur- 
bations do not necessarily inhibit fragmentation but are in- 
stead capable of inducing fragmentation in disc that otherwise 
would silently disperse or be accreted by the central star with- 
out ever experiencing fragmentation. The perturber star may 
accrete gas from the target circumstellar disc and may form 
planets on misaligned or even counter-rotating orbits with re- 
spect to the stellar spin. Future work will analyse discs with 
different masses and more or less massive host stars as well. 
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